Testing f{R) dark energy model with the large scale structure 
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In this work, we further investigate the family of ,f{R) dark energy models that can exactly 
mimic the same background expansion history as that of the ACDM model. Instead of using the 
parameterized framework of modified gravity, we study the large scale structure in the f{R) gravity 
using the full set cosmological perturbation equations. We investigate the structure formation in 
both the spatially flat and curved Universe. We also confront our model with the latest observations 
and conduct a Markov chain Monte Carlo analysis on the parameter space. 
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I. INTRODUCTION 

There has been accumulated conclusive evidence from 
supernovae[l] and other observations [2, 3] in the last 
decade, indicating that our Universe is undergoing a 
phase of accelerated expansion. Understanding the na- 
ture of the cosmic acceleration is one of the biggest ques- 
tions in modern physics. 

The leading explanation of the accelerated expansion is 
the cosmological constant within the context of General 
Relativity. However, the measured value of cosmolog- 
ical constant is far below the prediction of any sensible 
quantum field theories and the cosmological constant will 
inevitably lead to the coincidence problem that why the 
energy density of matter and vacuum are in the same or- 
der today (see [4] for review). Another possibility for the 
acceleration is that the universe is driven by a new and 
yet-unknown component called dark energy. The dark 
energy is some kind of dynamical fluid with negative and 
time-dependent equation of state w{a). However, the na- 
ture of the dynamical dark energy is even harder to be 
understood in the fundamental physics than that of the 
cosmological constant. 

Alternatively, a promising explanation of the accelera- 
tion is the modified gravity. The General relativity might 
not be ultimately correct in the cosmological scale. The 
Universe might be described by some kind of modified 
gravity. One simplest attempt is called f{R) gravity, in 
which the scalar curvature in the Lagrange density of 
Einstein's gravity is replaced by an arbitrary function of 
R{5]. The f{R) gravity can produce the accelerated ex- 
pansion of the Universe with any designed effective dark 
energy equation of state Furthermore, the time de- 
pendent effective DE EoS in the Jordan frame can be re- 
produced from the dilation of the inertial mass in the Ein- 
stein frame through the conformal transformation [7]. In 
this sense, the time dependent dark energy phenomenon 
can be better understood in a physical way in the frame- 
work of modified gravity. 

In this paper, we investigate a specific family of /(i?) 
models that can reproduce the same background expan- 
sion history as that of the AGDM model since it has been 
argued that a valid f{R) model should closely math the 
AGDM background [ti]. The family of /(i?) model con- 



tains only one more extra parameter than that of the 
AGDM model and furthermore, it does have the well- 
defined Lagrangian formalism in the spatially flat Uni- 
verse, which is valid for the whole expansion history of 
the Universe from the past to the future. The model is no 
longer simply a phenomenological model. Although this 
model has been studied by a number of work within the 
parameterized framework of modified gravity [9] [10], it 
still needs to solve the full set of cosmological perturba- 
tion equations to get more accurate results on the scale 
k dependent growth history of the Universe when con- 
fronted with upcoming high precision data in cosmologi- 
cal surveys. Therefore, in this work, instead of using the 
parameterized framework of modified gravity, we investi- 
gate the impact of f{R) models on the large scale struc- 
ture using the full set of covariant cosmological perturba- 
tion equations. We will confront our f{R) model with the 
latest observations and conduct a Markov chain Monte 
Garlo analysis on the parameter space. We will also ex- 
plore the spatially non-fiat case in the f{R) gravity, which 
has not been addressed in previous work [10] [11] [12]. 

This paper is organized as follows: In section II, we 
review the background dynamics of the Universe in the 
f{R) gravity. In section III, we present the scalar pertur- 
bation equations in the synchronous gauge for f{R) cos- 
mology and study the impact of f{R) gravity on the large 
scale structure formation. In section IV, we present the 
fittings results by confronting our model with the latest 
observations. In section V, we summarize and conclude 
this work. 



II. THE BACKGROUND DYNAMICS 

We work with the 4-dimensional action in the f{R) 
gravity [13] 

^ ^ i / ^'^^/^/(^) + / ^"^-^^"^ ' (1) 

where — SttG. We consider a homogeneous 
and isotropic background universe described by the 
Friedmann- Robertson- Walker (FRW) metric 
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where da^ is the conformal space-hke hypersurface with 
a constant curvature 
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The dynamics of the Universe in the f{R) gravity is de- 
scribed by[13] 
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where F = ^^^iS^, the dot denotes the time derivative 
with respect to the cosmic time t and p is the total en- 
ergy density of the matter which consists of the cold dark 
matter, baryon and radiation p = pc + Pb + Pr- P is the 
total pressure in the Universe. If we convert the deriva- 
tives in Eq.(4) from the cosmic time t to x ~\na , Eq.(4) 
can be recast into 



tF 



,ld\nE 
dx'^ ' ^2 dx 
_ dp 
~ iEdi ' 

where 



1) 



dF 
dx 



,d\nE 2K 



E 



-2X-- 



dx 



F 



(5) 



E ; 

R 

dp 
dx 



IP 

^0 



dx 



AE) 
3{p + p) . 



6Ke- 



(6) 



For convenience, the energy density p, K and the scalar 
curvature R in Eq.(5) are in the unit of Hq and we also 
set = 1 in our analysis. In order to get a viable f{R) 
model with a reasonable expansion history of the Uni- 
verse and without loss of generality, we can parameterize 
the quantity E(x) in Eq.(5) as the standard model in 
Einstein's gravity with an effective dark energy equation 
of state(EoS) u;[f)][14] 



where 
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= -K 
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indicates the energy density of the effective dark en- 



ergy, p 
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matter and p 
the matter. 



-p^ is the energy density of non-relativistic 
= Pj„ -|- p° is the total energy density of 



The f{R) model can be constructed by specifying the 
background expansion history of the Universe in Eq.(5). 
Eq.(5) becomes a second order differential equation only 
with respect to F{x). However, any designed background 
expansion history is meaningless in physics because we 
have less knowledge about the nature of the dark energy 
equation of stat w at present. Therefore, it is of great 
interest to investigate the f{R) model that can reproduce 
the same background expansion history as that of the 
ACDM model w = —1. Therefore in this work, we only 
focus on this family of f{R) models. We numerically find 
the solution of Eq.(5) in order to include the curved case 
51^ 7^ 0. We perform our numerical calculation starting 
from the deep matter dominated epoch(ai ^ 0.03). The 
curvature term K can be neglected at this epoch and the 
analytical solutions of Eq.(5) give rise to 
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dx 
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where p+ = ^^2^ ^™ ^® ^ dimensionlcss constant. 
In our work, we use the parameter Dm to characterize 
the family of f{R) models instead of using the parame- 
ter Bq proposed in [G] because, as presented in the ap- 
pendix VI B, Dm directly relates to the covariant param- 
eter vu in the Lagrangian formalism for this kind of /(i?) 
models in the spatially flat Universe. Since the f{R) 
model investigated in this work and ACDM model can 
only be distinguished in their perturbed space-time, in 
the next section, we turn to the cosmological perturba- 
tions theory. 



III. COSMOLOGICAL PERTURBATIONS 

The cosmological perturbation theory for the f{R) 
gravity has been well studied in [to]. [l^] has exten- 
sively presented the perturbation equations for a wide 
family of modifled gravities. In the appendix, we summa- 
rize the self-consistent covariant perturbation equations 
in the f{R) gravity. 

In this work, we use the synchronous gauge which 
is defined by -0 = and B = 0. The syn- 
chronous gauge is widely used in the Einstein-Boltzmann 
codes [f(i] [17] [18] [19] to calculate the temperature and 
polarization power spectra of the cosmic microwave back- 
ground anisotropy. 

The synchronous gauge is characterized by the follow- 
ing two parameters 
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where rjT refers to the conformal 3-space curvature per- 
turbation 
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The perturbed modified Einstein equations are given by 
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and 
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The scalar curvature perturbation SR is given by 
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With the perturbation equations, we perform our nu- 
merical analysis based on the public available Einstein- 
Boltzmann code CAMB [is]. In the appendix, we sum- 
marize the details on how to modify the code. In Fig.l, 
we show the full spectrum of temperature anisotropy 
Cf^ for a few representative values of £>,„. The cos- 
mological parameters used in the numerical process are 
taken from the WMAP 7-year best-fitted values for the 
ACDM model [2] n^K^ = 0.0227, l^^/i^ = 0.112, 
0.966, = 2.42 x lO'^f^A = 0.729, r = 0.085 . From 
Fig.l, we can see that the f{R) gravity only affects the 
lowest multipoles in the power spectrum while not change 
the acoustic peaks. This observation is consistent with 
what found in[(i]. The small / power is suppressed as the 
decreasing of the value of Dm- However, there is a turn- 
ing around point, roughly about Dj„ < — 0.37(Bo > 1-5), 
after which there is a prominent enhancement in the 
power as further decreasing the value of Dm- Numeri- 
cal results of our work agree with the already presented 
results within the parameterized framework of modified 
gravity [9] [10]. 

In Fig. 2, we present the linear cold dark matter power 
spectrum at redshift zero z = 0. The f{R) gravity 
changes not only the amplitude but also the shape of 
the power spectrum due to the scale k dependent growth 
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FIG. 1: The angular powerspectrum of the temperature cor- 
relation 



history [G]. However, the /(-R) model investigated in 
our work should not change the position of the Baryon 
Acoustic Oscillation(BAO) peak in the two-point corre- 
lation function in the real space because the sound hori- 
zon is only determined by the background cosmological 
parameters[20][21] and the family of /(i?) models inves- 
tigated in our work has the same background expansion. 
By doing the Fourier transformation 



?M = A / dkk'PL{k) 



sin(A;r) 
kr 



(19) 



in Fig 3, we present the two-point correlation function 
in the real space. We can sec clearly that although the 
shape of the correlation function has changed a lot for a 
few representative values of the position of the BAO 
peak does not change under the Fourier transformation, 
which shows that our numerical results are consistent 
with the theoretical prediction. 

In the spatially curved Universe, the curvature mainly 
shifts positions of acoustic peaks as shown in Fig. 4. Al- 
though in the lowest multipoles, the impact of spatial 
curvature has similar behaviors as that of the modified 
gravity £>„i, the high precision measurement of the acous- 
tic peaks in the CMB temperature angular powerspec- 
trum in combination with BAO and Hq[1] can put very 
tight constrains on the spatial curvature and the spatial 
curvature, therefore, could not have significant influence 
on the large scale structure in the f{R) gravity. 

Compared with the CMB temperature angular power 
spectrum C;"^^(see Fig.l), there are no significant imprint 
of f{R) gravity on the temperature polarization corre- 
lation spectrum and polarization auto-correlation spec- 
trum Cf^,Cf'^{scc Fig. 5). After the qualitative analy- 
sis of the impact of f{R) model on the large scale struc- 
ture, in the next section, we will present cosmological 
constrains on the f{R) model from latest observations. 
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FIG. 2: The linear cold dark matter powerspectrum. 
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FIG. 3: The two-point correlation function in the real space 



IV. CONSTRAINS FROM OBSERVATIONS 

The parameter space of our model is 

p = {nbh^,nrnh^,eA,nk,Hio''^A,in,,T,Dm) , (20) 

"where 9a is the angular size of the acoustic horizon, and 
As is amplitude of the primordial curvature perturbation. 
To avoid the instabilities in the high curvature region, 
it requires that Dm < [22]. Priors for cosmological 
parameters are listed in table I. For the CMB data, wc 
use the seven-year WMAP(WMAP 7) CMB temperature 
and polarization po"wer spectra [2] . We also take the mea- 
surement of the power spectrum of LRG from SDSS DR7 




FIG. 4: The impact of the curvature Slfe on the angular pow- 
erspectrum of the temperature correlation 




200 400 600 800 1000 1200 1400 1600 1800 2000 
t 



FIG. 5: The TE and EE angular powerspectrum. 



[23]. Ho'wever, we limit our considerations of the matter 
po"wer spectrum in a relatively linear scale and we only 
take the samples "within bins for k < 0.1/i/Mpc. We add 
the supernovac data [24] and the present day value of 
the Hubble constant from the measurement of Hubble 
Space Telescope(HST) Hq = 74.2± S.Gkms-^Mpc"^ [2.5] 
in order to put tighter constrains on the background cos- 
mological parameters. The global fitting results using 
above data sets are listed in Table II and HI. In the fiat 
Universe, "we find the constrains on the parameter Dm 
as \Dm\ < 0.709(95%CL) and \Dm\ < 0.711(95%CL) 
in the curved Universe. These results are equivalent 
to Bo < 3.86(95%CL) and Bq < 3.88(95%CL) if we 
use the parameter Bq. The result is consistent with 
[11]. Although the /(i?) gravity has significant impact 
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TABLE I: Priors for cosmological parameters 
0.005 < flth'' < 0.1 
0.01 < flch'^ < 0.99 
0.5 < e < 10 
-0.2 <Qk < 0.2 
0.01 < r < 0.8 
0.5 <ns < 1.5 
2.7 < ln[10^°As] < 4.0 
-1.2 < < 



on the shape of the matter power spectrum(see Fig. 2), 
the SDSS LRG matter spectrum ean not put very tight 
constrains on Dm because there are degeneracies be- 
tween Dm and other cosmological parameters [11]. When 
Dm ^ —0.6{Bo ^ 3), the CMB temperature angular 
power spectrum Cf"^ of the f{R) model goes back similar 
to that of ACDM model. Therefore, there arc two peaks 
in the likelihood distribution of Dm which are clearly 
shown in black lines in Fig. 6 and Fig. 7. Clearly, the 
combination of the data set CMB+SN+HST+MPK can 
not put very tight constrains on the f{R) models. We 
need to add additional data set. 

As pointed out in [11], the f{R) gravity can produce 
the anti-correlation in the Galaxy-ISW angular power 
spectrum. However, the measurement of the Galaxy- 
ISW correlation favors the positive correlation. There- 
fore, the Galaxy-ISW correlation data set can put very 
tight constrains on the f{R) model. We consider the 
cross correlation 

Cf = ^ J PdkPmik,0)W[>{k)W/{k) , (21) 

where the window functions W/ {k) and Wj^{k) are given 
by 

W/ik) = -To J dz^i'f - ^)ji[kx{z)] , 
W^f(fc) = J dzbiz)Il{z)D{z)ji[kx{z)] , (22) 

where x{z) is the looking back time, ji{x) is the spherical 
bessel function and Tq is the temperature of the CMB. 

In order to improve the performance in the numerical 
process, we use the Limber approximation 

- / k'dkji[kx]ji[kx'] « ^fc^ , (23) 
and Eq.(21) reduces to 

Cf^Tj,.>i^pJ^..) . (24) 

J X X 

where 

Pgiik,z)^^P^ik)Smik,z)^i^~^) , (25) 



and 

P^(fc) = ^(A)"=-i . (26) 

In this work, we adopt the Galaxy-ISW correlation 
data from [2(i] and use the public available ISWWLL 
code[2()] to calculate the likelihood. We have turned 
off the contribution of weak lensing in the original code. 
The results of the joint likelihood analysis with the data 
sets from CMB-hSN-hHST-hMPK-HglSW are shown in 
Table II and III. The marginalized ID and 2D likeli- 
hoods for interested parameters are shown in red lines in 
Fig. 6 and Fig. 7. The Galaxy-ISW correlation data set 
has improved significantly the constrains on the param- 
eter of Dm up to \Dm\ < 0.109(Bo < 0.376)(95%CL) in 
the flat Universe and \Dm\ < 0.131(Bo < 0.459)(95%CL) 
in the curved Universe. The result in the flat Universe are 
in good agreement with the work done by [].()]and [12]. 
The best fltted point for the curvature is slightly positive 
rife = 0.0063lo'Qog^ suggesting that the Universe is nearly 
flat but open {K < 0), which is diS'erent from the results 
that the Universe is nearly flat but closed(A' > 0) in the 
non-flat ACDM model Qk = ~0. 00231: a^ooltl-]- In the 
f{R) gravity, our Universe may have a completely differ- 
ent expansion scenario in the far future, if the Universe 
is not exactly spatially flat. 



V. CONCLUSIONS 

In this work, we have studied the impact of the spe- 
ciflc family of f{R) models that can reproduce the same 
background expansion history as that of ACDM model 
on the cosmic microwave background and the large scale 
structure using the full set of covariant cosmological per- 
turbation equations. Based upon the covariant pertur- 
bation equations, we have modified the public available 
Einstein-Boltzmann code CAMB and CosmoMC [27] to 
conduct the Markov chain Monte Carlo analysis on the 
parameter space and confront our f{R) model with the 
latest observations. In the flat Universe, our results are 
in a good agreement with previous independent work 
within the parameterized framework of modified grav- 
ity. We also have extended our analysis to the non-flat 
Universe. From the fitting results, we find that the Uni- 
verse is nearly fiat but open at the best-fitted point in 
the f{R) gravity which is different from the result that 
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TABLE III: Fitting results for the non-flat model fifc / 
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FIG. 6: Marginalized posterior distribution and 2-D contour plots for f{R) parameters in the flat Universe ^Ik 
CMB+SN+HST+MPK(black lines), CMB+SN+HST+MPK+gISW(red lines) 
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FIG. 7: Marginalized posterior distribution and 2-D contour plots for f{R) parameters in the non-flat Universe Q,k 7^ 0. 
CMB+SN+HST+MPK(black lines), CMB+SN+HST+MPK+gISW(red lines) 



the Universe is nearly flat but closed in the case of the 
ACDM model [2]. The modified gravity may change the 
final fate of our Universe in the future expansion, if the 
Universe is not exactly spatially flat. 
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curvature perturbations are characterized by Z and a 



h'L 
2k 

a ~ ka 



VI. APPENDIX 



A. Modifying CAMB 

where 

Our work is based on the public available Einstein- 
Boltzmann code CAMB [18]. CAMB propagates the 
covariant equations rest upon the CDM frame, which 
is equivalent to the synchronous gauge. Here we sum- 
marize the cosmological perturbation equations with the 

same conventions as that used in CAMB. In CAMB, the '7t 3 (f^ - 2^) • (27) 
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The perturbed modified Einstein equations can be writ- 
ten as 



^F{l3,a ~ Z) 
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where P2 is the curvature factor 
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' Fk F ■ 
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(28) 

(29) 
(30) 

,(31) 
(32) 



The propagation of the perturbed field 5F is given by 



5F" + 2'H5F' + a2(^ + m^)5F 
2 2 

'^{5p-i5p)-kF'Z , 



(33) 



and the perturbation of the scalar curvature 5R is given 
by 



a^5R = 2kZ' + QkZn - ^k^p2riT 



(34) 



We have replaced the original perturbation equations 
with the above set of equations in the original CAMB 
code. Another important part we need to modify is the 
source term of the CMB temperature anisotropy [28] [29] 
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where g ~ —ie ^ ~ aneaxE ^ is the visibility function 
and e is the optical depth. C is given by 



Q={\i2 + \e2) , 



(36) 



where /2, E2 indicate the quadrupole of the photon in- 
tensity and the E-like polarization respectively [29]. In 
the original code, a and Z are calculated from the per- 
turbation equations in the standard Einstein's gravity. 
In our work, a and Z are calculated from perturbation 
equations in the f{R) gravity. 



B. Lagrangian formalism 

In the spatially flat Universe filled with dust matter, 
the analytical solution of Eq.(5) reads 



F{x) - 1 = C^{e^-Y-2Fi[q-,p--,r_--e^^^] 
+ D„, (e3- )P+ 2 i^i [g+ , p+ ; r+ ; - e^- ^] 



(37) 



where the indexes are 

q~ 

r+ 
r_ 

P+ 



1 + V73 

12 
1 - 773 



12 



1 + 

1 - 
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6 
73 



6 

73 



12 

5 - V73 



12 



(38) 



The viable f{R) model should be in a "chameleon" type 
[30] [31] where in the past of the Universe, the model 
should go back to the standard model 



lim F{x) = 1 



(39) 



Cm, therefore, should vanish in the solution Cm ~ 0. 

The Lagrangian formalism for the corresponding /(i?) 
model is given by 



(35) /(i?) = i?-2A 



A 



i?-4A 



iF^ 



<i+,p+ - 1;''+;- 



A 



i?- 4A 



(40) 



where A, w are constant parameters. 

m = An(i?o - 4A)f+/(P+ - l)/AP+-i 



(41) 



When 07 = 0, A is just the cosmological constant. The 
above expression is a well-defined real function on the 
real axes in the physical range R > 4A, which is different 
from the results given by [32]. Moreover, the Lagrangian 
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formalism does not have the future singularity because 
f{R) is finite at the point of i? = 4A. 



lim /afi?) = 2A- 

fl-)-4A 

w4(-511 + 79y73)r(2/3)r(-r_) 
(-5 + ^/73)(-l + V73)(7 + ^/73)^(-p_)^(g+) 
« 2A- 1.256n7 



(42) 



where T is the Euler Gamma function. 

When R < 4A the expression of Eq. (40) becomes com- 
plex, which is clearly unphysical. 



O O J 

+ 4— V+|fcB' + ^fcHB-2-HV 
a 3 3 

+ 2ni!' - met)' - 24)" + a'(2^ + \e)] 

o 

+ F'[\kB + ij' + 2'Hi;~2(j)'] 
o 

a" 2 

+ + -k^ + 2K] 

a 3 



5F" - SF'n + 2ipF" 



= F[-k'^iP-k^ 



-k^E 



-E" 



^2knB + nE' 

kB'] 



(47) 



C. Cosmological perturbation in f(R) gravity 



F'{-E' + kB)-k^6F 



(48) 



The perturbed line element for the scalar perturbation 
can be written as [33] 

ds^ = a^[~{l + 2ijM''^)dT^ + 2BY}"UTdx' 

+ (1 + 2(t)Y^'^)-i,jdx'dx^ + EY,l;.'^dx'dx^] .(43) 



7ij in the spherical coordinate can be written as 








r2 

r'^sin?^ 



(44) 



--K^a^q = F[k(j)' - kUt}) + 2H'^ B 



a" 1, , 1 ,K, 
—B + -kE'- -E'—] 
a 6 2 k^ 

F'[nB- ^kijj] + ^kSF' 

IknsF - Ibf" , 

2 2 



(49) 



where q = {p+p)v. The above perturbation equations are 
self-consistent covariant equations. We can show that un- 
der the infinitesimal coordinate transformation, namely, 
Eq.(50), the perturbation equations could keep the same 
form. 



where we also include the spatially curved case. Pertur- 
bation equations in the spatially flat Universe can also 
be found in [34]. Y,Yj and Yij are the scalar harmonic 
function, which are defined by 



(A fc2)r(^) = , 



1y 



1 -r{s) 



(45) 



1 



fc2 ly ■ 3 

The perturbed modified Einstein equations arc given by 



sy/"^ = {b - ki° - p')yI^'^ 



-kl3-nf)Y^''> 
3 



EY, 



(50) 



(£; + 2fc/3)yA 

The perturbation for the scalar field 5F satisfies 

SF^JbrSR , 
, _ dF 

where the perturbation of scalar curvature 57? is given by 



knB + K{3(j}+^E) - Ik^E] 
2 b 

F'i'lB + 3m'-l4')-lnSF' 

, ,3 o 3 a" fc2 
2 2 a 2 



a'^SR ^ -12i' en-ip' + A(f>k^ - 2KE 

a 

+ 2^pk^ + imcj)' + \Ek'^ + 60" 
o 

- QBUk - 120A" - 2B'k . 



(51) 



Inserting Eq.(51) into Eq.(47) to eliminate 0" and us- 
(46) i]2g Eq.(46) to eliminate k'^4>, the equation governing the 
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behavior of SF gives rise to and 
6F" + 2nSF' + a2(^ + M^)SF 



3 i? = ^ + ^ . (54) 

+ F'(4HV'-30' + + , a? 

where 

^^' = J(7^-i?) , (53) 
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